Heralded quantum entanglement between two crystals 
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Quantum networks require the crucial ability to entangle quantum nodes [T]. A prominent example 
is the quantum repeater [2 4 which allows overcoming the distance barrier of direct transmission 
of single photons, provided remote quantum memories can be entangled in a heralded fashion. Here 
we report the observation of heralded entanglement between two ensembles of rare-earth-ions doped 
into separate crystals. A heralded single photon is sent through a 50/50 beamsplitter, creating a 
single-photon entangled state delocalized between two spatial modes. The quantum state of each 
mode is subsequently mapped onto a crystal, leading to an entangled state consisting of a single 
collective excitation delocalized between two crystals. This entanglement is revealed by mapping 
it back to optical modes and by estimating the concurrence of the retrieved light state [5]. Our 
results highlight the potential of rare-earth-ions doped crystals for entangled quantum nodes and 
bring quantum networks based on solid-state resources one step closer. 



Quantum entanglement challenges our intuition about 
physical reality. At the same time, it is an essential in- 
gredient in quantum communication [BJ, quantum pre- 
cision measurements [7] and quantum computing [5] [§] ■ 
In quantum communication, photons are naturally used 
as carriers of entanglement using either free-space or op- 
tical fiber transmission. Yet, even with ultra low-loss 
telecommunication optical fiber, the transmission prob- 
ability decreases exponentially with distance, limiting 
the achievable communication distance to a few hundred 
kilometers jTUj. A potential solution is to use quan- 
tum repeaters [5] based on linear optics and quantum 
memories [<2 [4] with which the entanglement distribu- 
tion time scales polynomially with the transmission dis- 
tance, provided entanglement between quantum memo- 
ries in remote locations can be heralded. In this con- 
text, one can more generally consider prospective quan- 
tum networks [T] where nodes generate, process and store 
quantum information, while photons transport quantum 
states from site to site and distribute entanglement over 
the entire network. An essential step towards the im- 
plementation of such potential technologies is to create 
entanglement between two quantum memories in an her- 
alded manner H] . 

Experimental observation of heralded entanglement 
between two independent atomic systems for quantum 
networks was achieved using cold gas ensembles involv- 
ing either two distinct ensembles [5j , or two spatial modes 
in the same ensemble [TTlfT5] (note that in ref. [15] the 
entanglement was postsclcctcd). Complete elementary 
links of a quantum repeater (based on heralded entangle- 
ment) have been implemented in cold gas systems [T61IT7] 
and with two trapped ions |18j . For the realization of 
scalable quantum repeaters, solid-state devices are tech- 
nologically appealing |19j . In this context, important 
results have already been obtained with rare-earth-ion 
(RE) doped crystals used as quantum memories. These 
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results include light storage times greater than one sec- 
ond [20] , storage efficiency of 69% [21] and quantum stor- 
age of 64 independent optical modes in one crystal [55] 
(see also ref. [13]). Recent achievements [2H[25] include 
the storage, in a single RE-doped crystal, of photonic 
time-bin entanglement generated through spontaneous 
parametric downconversion (SPDC). Heralded entangle- 
ment between two RE-doped crystals using an SPDC 
source, a common resource in quantum optics, repre- 
sents an important step towards the implementation of 
quantum repeater architectures based on solid-state de- 
vices [55] . 

Here we present the observation of heralded quantum 
entanglement between two neodymium ensembles doped 
in two yttrium ortho-silicate crystals (Nd 3+ :Y2SiOs) sep- 
arated by 1.3 centimeters. Let us conceptually detail 
our experiment, depicted in Fig. 1. A nonlinear optical 
waveguide is pumped to produce photon pairs through 
SPDC. The resulting idler (1338 nm) and signal (883 nm) 
photons are strongly filtered to match the absorption 
bandwidth of the crystals, yielding a coherence time of 
7 ns. In the limit where the probability of creating a sin- 
gle pair is much smaller than one, the detection of an idler 
photon heralds the presence of a single signal photon. 
By sending the latter through a balanced beamsplitter, 
one ideally heralds a single-photon entangled state [28 
^(|1)a|0)b + |0)a|1)b) between the two spatial output 
modes A and B. In each of the modes a crystal acts as 
quantum memory, denoted Ma or Mb ■ Upon absorption 
of the single photon by one of the crystals the de- 
tection of an idler photon heralds the creation of a single 
collective excitation delocalized between the two crystals, 
ideally written as ^(\W) A \0) B + \0) a\W) B ). Here, \W) A 
(or |W)s) is the Dicke-like state created by the absorp- 
tion of a single photon in Ma (or Mb). To determine the 
presence of entanglement between the memories, we use 
a photon echo technique based on an atomic frequency 
comb [2"9"H3"3"] that reconverts the collective excitation into 
optical modes A and B after a preprogrammed storage 
time of 33 ns. Note that the memories are used in a 
double-pass configuration, which effectively increases the 
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FIG. 1. Experimental setup. The quantum memories Ma and Mb are implemented using neodymium ions doped into 
yttrium ortho-silicate crystals (Nd 3+ :Y2SiOs) separated by 1.3 cm and cooled to 3 K using a cryostat (see ref. |24j for details). 
The total efficiency of each memory was 15%. A fiber optic switch is used to alternate between a 15 ms long preparation of the 
two Nd ensembles as atomic frequency combs on the 4 F 3 / 2 — > 4 Ig/2 transition, followed by attempts of entanglement creation 
for another 15 ms. For the latter, continuous wave light at 532 nm is coupled into a periodically poled KTP waveguide, leading 
to the production of pairs of photons at wavelengths of 883 nm and 1338 nm through spontaneous parametric downconversion. 
Photons from each pair are separated on a dichroic mirror (DM) and frequency filtered to below the 120 MHz bandwidth of the 
quantum memories. The detection of an idler photon at 1338 nm (using a low-noise superconducting single photon detector [27]) 
heralds the presence of a signal photon at 883 nm. The signal photon now traverses the switch, a polarizing beam splitter (PBS) 
and a Faraday rotator (FR), before a 50/50 beamsplitter (BS) creates single-photon entanglement between spatial modes A 
and B. This entanglement is, upon absorption, mapped onto the crystals Ma and Mb- After a preprogrammed storage time 
of 33 ns, the photons are reemitted and pass through the BS again. Depending on which output mode of BS they emerge from, 
they either reach detector 1 or are rotated in polarization by the FR and reflected by the PBS towards detector 2. 



optical depth and the quantum memory efficiency up to 
15%. The resulting fields can then be probed using sin- 
gle photon detectors to reveal heralded entanglement be- 
tween the memories. Since the entanglement cannot in- 
crease through local operations on the optical modes A 
and B, the entanglement of the retrieved light fields pro- 
vides a lower bound for the entanglement between the 
two memories. 

The photonic state retrieved from the memories is de- 
scribed by a density matrix p, including loss and noise, 
expressed in the Fock state basis. To reveal entangle- 
ment in this basis, one cannot resort to violating a Bell 
inequality given solely inefficient and noisy single pho- 
ton detectors. Instead, as shown in ref. [3] (see also 
ref. [Ml E]), a tomographic approach based on single 
photon detectors can be used. Indeed, it is possible to 
determine the presence of entanglement between the re- 
trieved fields from the knowledge of the heralded prob- 
abilities Prrm of detecting m photons in mode A and n 
in mode B, where m,n € {0, 1}, and from the coherence 
between these modes. More precisely, one can obtain a 
lower bound on the concurrence C of the detected fields, 
a measure of entanglement ranging from for a separable 
state to 1 for a maximally entangled state, through 

C > max(0, V(p i + Pio) - VPooPn)- (1) 

The term V is the interference visibility obtained by re- 
combining optical modes A and B on a 50/50 beamsplit- 
ter: it is directly proportional to the coherence between 
the retrieved fields in modes A and B. To obtain a large 
concurrence, that is, a large amount of entanglement, one 
should maximize V (the coherence) and pio + Poi (the 



probability to detect the heralded photon), and mini- 
mize poo and Pn (the probabilities of detecting separable 
states 10)^10)5 and |1)a|1)b stemming from a lost sig- 
nal photon and from two signal photons, respectively). 
Implicit to this method is the assumptions that i) the 
creation of more than two pairs is negligible, and ii) the 
off-diagonal elements of p with different number of pho- 
tons vanish (this is valid since no local oscillator provid- 
ing a phase reference was used; see ref. [S]). 

To estimate V, poo, pio, poi and pn, we used the setup 
of Fig. 1 in the following way. First, the visibility V 
was measured by allowing the re-emitted delocalized pho- 
ton to interfere with itself using a balanced Michelson 
interferometer whose phase was actively stabilized (see 
Fig. 1). Then, we blocked spatial mode B to estimate 
Pio by summing the number of detections on detectors 1 
and 2, conditioned on a heralding signal. Probability poi 
is estimated similarly by blocking mode A. Note that 
the probability of getting simultaneous counts on both 
detectors was negligibly small in front of the probability 
of a single detection. Then, poo is estimated through nor- 
malization of the total probability, poo + pio + Pai ~ 1, 
which is justified by the fact that pn <C Pio + Poi Poo 
(see measured values in the Appendix). To estimate pn, 
we used two different methods, described below. 

In the first method, we use a direct measurement of 
threefold coincidences, i.e. involving all three detectors. 
Details about this measurement are given in the Ap- 
pendix. With a pump power of 16 mW and a coin- 
cidence time window of 10 ns, we obtained C^ MLE ) = 
6.3(3.8) x 10 -5 using a maximum likelihood estima- 
tion (MLE) of the threefold coincidence probability, and 
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(7( CE ) = 3.9(3.8) x 1CP 5 using a more conservative es- 
timation (CE); see Appendix. Both estimations yield a 
concurrence that is greater than by at least one stan- 
dard deviation, and show that entanglement was indeed 
present between the two crystals. This measurement re- 
quired 166 hours, a period in which two threefold coin- 
cidences were observed. The prohibitively long integra- 
tion time of this method prevented us from attempting 
it with lower pump powers (i.e. for lower probability of 
creating more than one pair). Hence, to study how the 
concurrence changes with pump power, we used a second 
method based on twofold coincidences, which we now de- 
scribe. 

In the second method, pu is estimated using a sup- 
plementary assumption (see the Appendix for details). 
This approach is motivated by the results obtained in 
ref. [13]. Specifically, we assume that all the observed 
detections stem from a two-mode squeezed-state, and 
thus the measured zero-time cross-correlation g s ^ can be 
written as g s ,i = 1 + 1/p, where p « 1 and p 2 are in- 
terpreted as the probabilities of creating one and two 
photon pairs, respectively. We then proceed as follows 
to estimate p\\. We first measure the zero-time cross- 
correlation gf^ i between detections in the idler mode and 
the signal mode A (with mode B blocked). Then we 
measure gf i in the same way (with mode A blocked) and 
verify that gf i « gf i . We calculate the average of gf i 
and gf i} denoted g s ^, and estimate pu using 



In the Appendix, we provide justifications and additional 
measurements that support our assumption and give ev- 
idence that it yields a lower bound on the concurrence. 
In particular, we measured the second-order autocorre- 
lation of the signal (or idler) mode without storage to 
be 5 i 2 j(0) = 1.81(2) (or g$(0) = 1.86(9)), which is very 
close to the ideal value of 2 associated with the ther- 
mal photon statistics of a two-mode squeezed state. We 
also measured the zero-time second-order autocorrelation 
function of the heralded signal photon just before stor- 

(2) 

age and obtained g% {i (0) = 0.061(4) for a pump power 
of 8 mW, which is consistent with p <C 1. Note that 
this method of estimating pu is based on twofold coinci- 
dences as opposed to threefold coincidences for the other 
method. For our specific setup, this resulted in a reduc- 
tion of the measurement time by a factor of 10 6 for similar 
statistical confidence on the concurrence. This method 
requires no physical modifications to the optical circuit to 
measure the different components of the retrieved fields, 
which simplifies its implementation. 

We performed a series of measurements for several val- 
ues of the pump power, which is proportional to p pro- 
vided p -C 1. For all measurements we used coincidence 
windows of 10 ns, and all results are based on raw counts 
(i.e. without subtraction of dark counts and accidental 
coincidences). The inset of 2a shows the interference of 
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FIG. 2. Results, a, Visibility as a function of pump power. 
The visibility is approximately constant with an average of 
96.5 ± 1.2% (green shaded region). The inset shows the visi- 
bility curves at 16 mW measured with detectors 1 (•) and 2 
(□) (15 minutes acquisition time per point; the error bars are 
smaller than the symbols). The different amplitudes result 
from non-uniform loss after recombination of modes A and B 
on the beamsplitter. The visibilities agree within uncertain- 
ties in the fits, b, Zero-time cross-correlation g St i as a function 
of pump power. The decreasing values agree well with a theo- 
retical model (shaded region; see Appendix), c, Lower bound 
on the concurrence estimated using the cross-correlation mea- 
surement the as a function of pump power (•) , calculated us- 
ing Eq.[I]and[2] and the mean visibility of a. The concurrence 
decreases with pump power, as expected, but remains positive 
up to 16 mW. The shaded region corresponds to the model 
in b. All measured values of pio, poi and pu are given in the 
Appendix. All values are based on raw counts. Uncertainties 
are obtained assuming Poissonnian detection statistics. The 
lower bound on the concurrence at 16 mW obtained from 
measured threefold coincidences using either the maximum 
likelihood estimation (O) or the conservative estimation (♦) 
are also shown (they are horizontally offsetted for clarity). 
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the delocalized single photon retrieved from the memo- 
ries as a function of the phase of the interferometer. The 
visibilities obtained with detectors 1 and 2 agree within 
the uncertainty of the fit. As shown on Fig. 2a, the vis- 
ibility does not depend on the pump power and has an 
average of 96.5(1.2)%, implying almost perfect coherence 
between the two memories and a large signal-to-noise 
ratio. Fig. 2b shows the measured g s .i- Reducing the 
pump power increases the cross-correlation, as expected. 
Using additional measurements, we estimated the trans- 
mission loss, memory efficiency, dark count probability 
and the pair creation probability of our setup. These val- 
ues were then used in a theoretical model (shaded region 
on Fig. 2b) that is in excellent agreement with the mea- 
sured values of g S i. As explained in the Appendix, this 
agreement provides evidence that the estimated value of 
Pn yields a lower bound on the concurrence. Finally, 
Fig. 2c shows the lower bound on the concurrence for 
all pump powers, calculated using Eqs. [T] and [2] and the 
average visibility of Fig. 2a. The concurrence decreases 
with pump power because pn increases, while all other 
terms of Eq. [T] depend on photon loss only and hence are 
approximately constant. Nevertheless, the concurrence 
stays positive for all used pump powers, proving the cre- 
ation of heralded entanglement between the atomic en- 
sembles inside the two crystals. The theoretical curve on 
Fig. 2c (shaded region) is based on the same model used 
for the theoretical curve in Fig. 2b. 

The results of the measurement of the concurrence 
based on threefold coincidences are plotted on Fig. 2c and 
agree, within uncertainty, with the result of the method 
based on measurement of the cross-correlation. 

The observed values of the concurrence are lower 
bounds on the amount of entanglement of the detected 
fields and are almost entirely determined by optical loss. 
To see this, we first note that when both the multi-pair 
creation probability and the total transmission probabil- 
ity i] of the signal photon are small, the concurrence is 
approximately given by C ~ rj(V — 2/ yJg St i — 1). At 
8 mW of pump power, we measured g s .i ~ 10, poo = 
0.9997831(71), and p u = 5.18(40) x 10~ 9 , and hence 
C fa 0.3?7. With 77 « pio + poi = 2.2 x 10" 4 , we see 
that C R3 6.6 x 10~ 4 <C 0.3. Using this, an estimate of 
the concurrence of the fields retrieved just after the crys- 
tals is obtained by noticing that the transmission r\ is 
the product of the probability to find the signal photon 
inside the fiber per heralding signal (20%), the memory 
efficiency (15%), the transmission in the interferometer 
(2.4%) and detector efficiency (30%). Subtracting de- 
tector inefficiency and interferometer loss yields a lower 
bound of approximately 0.092 for the concurrence of the 
field retrieved just after crystals. 

The high interference visibility indicates that the 
stored entanglement does not significantly decohere dur- 
ing the 33 ns, and the coherent nature of the storage 
endures for longer times, as shown previously using stor- 
age of weak coherent states [22J [3TJ] . Increasing the stor- 
age time does however lower the cross-correlation func- 



tion g s ,i, as measured in ref. [23], and this should conse- 
quently lower the concurrence of the stored entanglement 
(this is essentially due to the decrease of the memory effi- 
ciency with increasing storage times, as shown in the Ap- 
pendix and in ref. [23]). With Nd 3+ :Y 2 Si0 5 , one could in 
principle store entanglement up to ~ 1 us (see ref. [22]). 
A promising approach to go beyond these limits, and to 
allow on-demand retrieval of the stored entanglement, is 
to implement spin-wave storage [29] to increase storage 
times (as demonstrated recently [31] in Pr 3+ :Y 2 Si05), 
and to place the crystal inside an impedance-matched 
cavity to increase the efficiency |34j . Such improvements 
are necessary for the development of a quantum repeater 
based on solid-state quantum memories, and are the sub- 
ject of current research. 

In conclusion, we have reported an experimental ob- 
servation of heralded quantum entanglement between 
two separate solid-state quantum memories. We em- 
phasize that although the entangled state involves only 
one excitation, the observed entanglement shows that the 
stored excitation is coherently delocalized among all the 
neodymium ions in resonance with the photon, meaning 
roughly 10 10 ions in each crystal. Our results demon- 
strate that rare-earth-ion ensembles, naturally trapped 
in crystals, have the potential of compact, stable and co- 
herent quantum network nodes. Moreover, single photon 
entanglement is a simple form of entanglement which can 
be used for teleportation [35] and entanglement swapping 
operations [3] and can also be purified using linear opti- 
cal elements [33]. It is also a critical resource in several 
proposals for quantum repeaters [3] [26l [37] since it is 
less sensitive to transmission loss and detector inefficien- 
cies [3] as compared to two-photon entangled states. 

Our experimental approach is based on solid-state de- 
vices, the key components being the PPKTP chip with an 
integrated waveguide (the photon source) and the crystal 
memory. We believe this approach opens possibilities of 
integration of components, for instance frequency filter- 
ing directly on the chip [55] , or waveguide quantum mem- 
ories [25]- The prospect of combining solid-state photon 
sources and quantum memories is attractive for practical 
future quantum networks. One important challenge in 
this context is to create entanglement between two re- 
mote solids in a heralded way, which would realize an 
elementary link for quantum repeaters. 
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APPENDIX 

This appendix provides details on the method we used 
to estimate the threefold coincidence probability pn us- 
ing either the cross-correlation g s ,i or by direct measure- 
ment. Specifically, in section |Aj we present a model de- 
scribing our setup composed of the source of photon pairs 
and the memories, and includes the effect of dark counts. 
We also detail additional measurements that support our 
model. We present evidence that the method based on 
the measurement of g Si i to estimate pn yields a lower 
bound on the concurrence at detection and, consequently, 
on the concurrence of the delocalized excitation stored in 
the quantum memories. In section [B] we provide details 
on our method to estimate p\\ directly, we show how 
the uncertainty on the threefold coincidence probability 
is calculated and we show the results. 



Appendix A: Concurrence using the 
cross-correlation 



Let us begin by a description of the source. In our 
model, the SPDC source produces a two-mode squeezed 
state (where the coherences are lost due to the lack of a 
phase reference) 



+00 

cosh -2 r tanh 2 ™ r \nin 8 ) {riin s 
n=0 



(Al) 



Here, and n s correspond to the number of photons 
in the idler and signal modes (and are both equal to n). 
While the signal and idler fields are perfectly correlated in 
photon number, they individually exhibit thermal statis- 
tics and hence, their second-order auto-correlation func- 
tion g( 2 \0) is, in theory, equal to 2. We experimentally 
tested this for the photons emitted by the source by 

(2) (2) 

performing direct measurements of (t), gsj(r) and 

The latter is the auto-correlation 
conditioned on the detection of an 



A3 



function of the signa 



idler photon at time r. For the measurements we added 
50/50 beam splitters after the frequency filters of the 
signal or idler photon, respectively, and registered coin- 
cidence events using free running detectors. We found 

(2) 

zero-delay auto-correlations of 9^/(0) = 1.86(9) and 

(2) 

9sJ(0) = 1.81(2), which are very close to the expected 
values. This is consistent with the assumption that our 



source produces a state of the form of Eq. (All. Addi 



tionally, we found a conditional auto-correlation of the 
signal photon at zero time delay of g s ^(0) = 0.061(4) at 
a pump power of 8 mW. This shows the very good single 
photon character of the heralded signal fields, and that 
tanh 2 r < 1. 

Let us continue the description of our model by fo- 
cusing on the memory. Due to non-unit absorption effi- 
ciency, each signal photon detection after the memory 
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FIG. A3. Measurements of second-order auto-correlation 
bunctions of a, the idler photon (1338 nm), b, the signal 
photon (883 nm) before it is sent towards the crystals, and c, 
the same as b but conditioned on a coincidence detection of 
an idler and the first signal photon at delay r = 0. The values 
at zero delay in a and b are close to the theoretical value of 
2. The almost zero minimum value in c gives evidence of the 
single photon character of the heralded signal photon. Data 
was accumulated for 72 hours at a pump power of 16 mW in 
a, and for 12 hours at 8 mW in b and c. The values in c have 
been calculated for a coincidence window of 10 ns. Error bars 
are due to Poisson counting statistics. 



stems either from a photon created at an early time, 
stored in the medium, and later retrieved; or from a 
photon created at a later time and not absorbed by the 
medium. Hence, we need to consider two time bins. The 
storage process is a linear operation, mapping the light 
field onto an atomic state. Thus, in our model, it can be 
represented by a beamsplitter with a transmission cor- 
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FIG. A4. Model of our experimental setup. The source pro- 
duces idler-signal pairs described by a two-mode squeezed 
state. The memory is represented by three beamsplitters 
(to account for the absorption, the decoherence and the re- 
emission respectively) . Each temporal mode after the memory 
d a comes either from a signal photon emitted at an early time, 
then stored and later retrieved (mode s e ), or from a photon 
produced at a late time and directly transmitted through the 
memory (mode s^). To take both contributions into account, 
a virtual source creating idler-signal pairs in a late time bin 
has been introduced. The correlations between the modes 
d s and di are characterized through the measurement of the 
second-order cross-correlation function. 



responding to the absorption probability, c.f. Fig. |A4| 
Another beamsplitter is introduced to account for the 
decoherence, inherent to the storage process [53]. Fi- 
nally, a last beamsplitter accounts for the reemission of 
the photon. It is identical to the first one and combines 
the two time bins into a single temporal mode d s . The 
output of the last beamsplitter is detected to access the 
cross-correlation function 



9 s 



(dldi4d s 



(d]d z )(dld 



(A2) 



between detection on the early idler mode di and signal 
mode d s . Note that the overall transmission (in intensity) 
of the three beamsplitters corresponds to the efficiency 
of the storage and retrieval processes ?7ccho, while the re- 
flection of the last beamsplitter is the transmission of the 
atomic ensemble ?7trans- 

Using the Heisenberg picture 



where 



(d\d s ) = (0\UUlUUU s U\0) 

Wd s U= V^ocho (cosh r s e + sinhr i\) 
+\/ 7 7traiis (coshr s e + sinhr if] 

leading to 

(d\d 3 ) = (rjtrans + Vecho) sinh 2 r. 
Furthermore, we have 

U'diU — (coshr i e + sinhr s\) 

leading to 

(d\di) = sinh 2 r. 



Following similar lines, we find 

{d\did\d s ) = sinh 2 r( sinh 2 r (rfeans + r? echo ) 

+ COsh 2 r 7/echo) • 

Therefore, the cross-correlation function between the 
idler-signal modes is given by 



1 + 



f 



tanh r ( 1 



Note that the effect of the memories is to add the term 
(1 + atiasa), which decreases the cross-correlation. One 
can further take the detector noise into account by adding 
an additional source of noise with a Poissonian photon 
distribution. We finally find 



1 + 



tanh r 



(A3) 



where ?7dark is the probability to get a dark count within 
the detection window and p c is the conditional proba- 
bility of detecting the signal photon retrieved from the 
memory. Note that for small pump powers, the proba- 
bility to create one pair is given by 



tanh 2 rj cosh 2 



OiP 



pump ) 



where P pU mp is the continuous wave pump power. 

To support the model presented above, we performed 
additional measurements to determine the values of a, 
?7trans/f?echo) Pc and (see Table Al). First, we estimated 
a = 2.71(8) x 10~ 3 pairs/mW from a measurement of 
the cross-correlation as function of pump power with 
one crystal prepared as a transmission window (no stor- 
age) while the path of other crystal was blocked. Then, 
^trans/^ccho was measured from the ratio of the proba- 
bilities of detecting the heralded signal photon when it 
is either transmitted through the memory (no storage), 
or stored and retrieved from the memory. This measure- 
ment was performed twice - once for each arm of the 
interferometer (with the other arm blocked) - and the 
outcomes were averaged. The conditional probability p c 
was taken as the average of probabilities p w and poi of 
detecting the heralded signal photon retrieved from the 
memories (see main text). Finally, we measured a dark 
count probability of 7?dark = 2 x 10~ 6 per 10 ns detec- 
tion window for the signal mode detector. Idler photons 
were detected using a superconducting nanowire single 
photon detector with a dark count rate of 20 Hz. This 
low dark count rate is neglected in our calculations. Us- 
ing these parameters, we compared the measured cross- 
correlation g s i with values of g Sj i predicted by Eq. (A3 1. 



The comparison is shown on Fig. 2B in the main text. 
We find excellent agreement within the statistical uncer- 
tainties, thereby supporting our model for the source and 
the memory. 
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TABLE Al. Conditional probabilities poi, Pio, Pn and the 
ratio ?7trans/?7ccho for various pump powers. 



Pump 


poi (1(T 4 ) 


pio (10 -4 ) 


Pn (io -9 ) 


%rans 


power 


lecho 


1 mW 


1.04(14) 


0.82(12) 


1.33(30) 


2.84(33) 


2 mW 


1.193(75) 


0.809(63) 


1.63(19) 


3.03(17) 


3 mW 


0.952(72) 


0.878(70) 


1.61(20) 


2.59(17) 


4 mW 


1.105(72) 


0.902(66) 


2.82(31) 


3.35(19) 


8 mW 


1.185(51) 


0.984(50) 


5.18(40) 


3.13(12) 


13 mW 


1.247(56) 


1.131(52) 


8.79(66) 


2.86(11) 


16 mW 


1.146(47) 


1.175(48) 


9.56(64) 


2.748(93) 


Mean 


1.123(30) 


0.957(27) 




2.936(69) 



We now give evidence based on our model that the 
method we used to estimate the threefold coincidence 
probability leads to an underestimation of the concur- 
rence. For this, we compare the threefold coincidence 
probability i predicted by our model with the one ob- 
tained from the method based on the measured cross- 
correlation g s j. Let us first concentrate on p*^. Using 
our model, we can predict that the threefold coincidence 
probability should be given by 



Pn = 4pioPoi 



OtP pump I 1 



?7tra 



echo 



?7dark 
Pc 



This should be compared with the threefold coincidence 



probability estimated using g s 



1 + 1/p, i.e. assum- 



ing detections stem from a two-mode squeezed states 
(see mai n text). Combining this method with our model 
(Eq. [A3} yields 

%>loffoi 

'/-.. - 1 



Pn 



4pioPoi 



aP, 



pump 



l 



7/ccho 



'Mark 
Pc 



One can see that the ratio Tftrans/^echo enters inpn with- 
out the prefactor of one-half, such that p u > From 
Eq. (1) of the main text, we directly see that an overesti- 
mation of pn leads to a lower bound on the concurrence 
at detection, and hence on the concurrence of the entan- 
glement stored in the memories. 



Appendix B: Concurrence using a direct 
measurement of the threefold coincidence 
probability 

1. Method 

We now describe how we estimate the probability pn 
from the measurement of threefold coincidences. For this, 



we use the experimental setup shown on Fig. A5 this is 
the same setup shown on Fig. 1 of the main text. Unlike 
our estimation based on the measurement of the cross- 
correlation (section [A]), this direct method does not re- 



quire assumptions on the statistics of the source of pho- 
ton pairs and on how the memories operate. 

Let us consider the two-photon part, denoted p 2 , of the 
complete density matrix describing the state of the re- 
trieved fields just after the crystals. If the relative phase 
between the spatial modes is randomized, then the co- 
herences are zero and we can write 

P2 = <Zn|ll}(ll| + <Z2o|20)(20| + g 02 |02)(02|. (Al) 

Here, qu is the probability per heralding signal that the 
fields retrieved from crystals A and B each contain one 
photon, (720 (902) the probability per heralding signal that 
the fields retrieved from crystals A and B each contain 
2 and (0 and 2) photons, respectively. Estimation of 
Pn ca n be performed by measuring the 911 probability 
in Eq. Al as in Ref. [5]. We show here that one can also 
estimate pn by measuring g 2 o an d 902- 

We first derive useful relations between the occupa- 
tion probabilities. We use the following definitions: Pjj 
is the probability to herald the presence of a signal pho- 
ton (i.e. the probability to detect an idler photon dur- 
ing the detection time window); rjf cho (Vecho) * s * ne stor- 
age and retrieval efficiency of crystal A (B); P2 is the 
probability that the source of photon pair emits 2 pairs 
of photons during the detection time window; R and T 
are the intensity reflection and transmission coefficients 
of the beamsplitter when photons are incident from the 
left hand side, satisfying R + T < 1 and R, T > (see 
Fig. A5). Using these definitions, and the linear nature 



of the optical storage process, we have 

911 = TT 2RTl 1echoVecho P 2, 

920 = -^T 2 ( V i ho ) 2 P 2l 

920 - ^R 2 {VechofP2, 

with which we can write 920 and 9q 2 as a function of qu: 



920 
902 



9ii(i?/2T), 
qu(T/2R). 



(A2) 



We use these expressions later. 

We now consider recombination at the beamsplitter. 
We describe the effect of the beamsplitter on mode op- 
erators a and b as follows: 

& PrO* + fob* + , 



where mode I represents loss. We set all coefficients to be 
real. Energy conservation yields a 2 + a 2 + of = 1 (and 
similarly for the coefficients of mode b). Also, unitarity 
imposes that a r at + firfit = 0. Finally, we suppose that 
< a 2 ,a 2 ,f3 2 ,/3 2 < 1/2, and we note that R — f3 2 and 
T = a 2 . Using these notations, we can show that the 
probability to get a coincidence at detectors 1 and 2 per 
heralding signal is given by 



P11 = (an9n + «2o920 + 002902)7/1 f?2, 



(A3) 




FIG. A5. Experimental setup used to estimate pn directly. This setup is the same of Fig. 1 of the main text. The relative 
phase between the arms of the interferometer was randomized by letting the interferometer (see Fig. A5 I drift by itself over the 
total measurement period of more than 100 hours. PBS: polarizing beamsplitter; FR: Faraday rotator; BS: beamsplitter; Ma 
and Mb are crystal memories A and B. The detection efficiency from the output modes of the beamsplitter to detectors 1 and 2 
and r\\ and 772, respectively. The amplitude transmission and reflection coefficients of modes a and b right after the crystals, 
through the beamsplitter, are shown. 



where 771 (772) is the transmission from output mode a 
(b) - after the beamsplitter - to detector 1 (detector 2), 
including the detector efficiency. We also have 



an = Pr 



a 2 o = 2a^al, 
a 02 = 2/3ftf. 



01 



We characterized our beamsplitter and obtained oq 



T = 0.479, a z r = 0.422, 
These values yield 



0.482 and # = R = 0.409. 



an = 0.0028, 
a 20 = 0.394, 
a 02 = 0.404. 



(A4) 



We immediately notice that an can be neglected with re- 
spect to <i2o and 002- This essentially means that when- 
ever a photon is retrieved from each crystal, they bunch 
with almost certainty and do not significantly contribute 
to threefold coincidences. It also requires indistinguisha- 
bility of the retrieved photons, which is consistent with 
the large single photon interference visibility observed 
(96.5%) 



O20 



Combining Eq. 
a 2 ~ 0.4, we get 



A3 and Eq. A4 and writing 



Q20 + 102 



pn w 0.4 7/1 772(1720 + 902) 

Finally, using Eq. |A2| we get 

R T 
2T + 2R 
= 1.012 q n 

~ qn 

and thus 

pn bs 0.4771772911. 



(A5) 



9ii 



(A6) 



We now have to determine how to use Eq. |A6| in the 
calculation of the concurrence. This method should give 



a result consistent with the way we measured p±o + P01, 
which was done as follows. We let the phase of the inter- 
ferometer drift to randomize it, and measured the prob- 
ability to get a detection at detector 1 or 2 per heralding 
signal. This is equivalent to, first, measure crystal A 
(with the path of crystal B blocked) using an effective 
detector of efficiency 77^ = 0^771 + 0:^772, and then mea- 
sure crystal B (with the path of crystal A blocked) using 
an effective detector of efficiency 77s = PrVi + Pt^i and 
then add the results together. Hence, if both crystals 
were measured simultaneously and separately, the three- 
fold coincidence probability per heralding signal would 
be 



(A7) 



qnVAVB- 



We can now find the relation between Eq. |A6| and |A7| 
For this, we use the fact that, in our experiment, 772 was 
larger than 771, which is apparent from the visibility plot 
in the main text (inset of Fig. 2- A). Hence, by setting 
771 = 772/2, we get qnr]Ar]B ~ 0.454 rj^qii- Therefore, we 
immediately see that if we multiply Eq. |A6| by 2.27, we 
get Eq. \K7\ 

In summary, to estimate the probability pn used in 
the formula of the concurrence, namely 



C > max(0,V^(pio + P01) - 2VPooPn), 



(A8) 



we multiply the measured pn by 2.27. This multiplica- 
tive factor was obtained using approximations, and hence 
has some uncertainty associated to it. Nevertheless, as 
explained in the next subsection, this uncertainty can 
be neglected in front of the large statistical uncertainty 
of pn. 



2. Statistical uncertainty 

Here we provide details on how the statistical uncer- 
tainty on pn was calculated in the method based on the 
measurement of threefold coincidences (subsection Bll. 
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We denote Nh the number of heralding signals during 
the measurement period T, Pn the stochastic variable 
representing the probability to register a threefold coin- 
cidence per heralding signal, and pn our estimation of 
Pn obtained from the measurement. 

Since pn <C 1 and Nh ^> 1, the distribution of the 
number of coincidences between detectors 1 and 2 per 
heralding signal (i.e. threefold coincidences) is accurately 
described by a Poisson distribution. Specifically, the 
probability to get n coincidences given Nh heralding sig- 
nals is given by 



P{n\v\i) = e 



- NhPii (NhVuT 



(A9) 



The probability P(n\pn) is the likelihood function of pn 
given n. Hence, Pn can be estimated by maximizing this 
likelihood, which yields 



(MLE) _ n yn 
Pn ~N H ± N H ' 



where MLE stands for "maximum likelihood estimation" . 
This method yields an appropriate estimate of Pn , unless 
one measures n = 0, in which case one would conclude 
that p n = ± 0. 

A more conservative estimate can be obtained in the 
following way. Using Bayes' law, we can write the prob- 
ability density function fl(pn\n) of getting Pn — pn 
given n measured threefold coincidences: 



0(pn|n) 



P(n\pi 



Jo P(n\Pu)dpu 
The denominator of Eq. |A9| can be calculated 

1 



-N. 



\- N ^ iNH ^dp 11 ^ — ll-e _ „ 
n\ N H \ ^ fc! 

1 

~ N~h~' 



The last equality holds because n <C Nh ■ Therefore: 

fi(pii|n) = N B P(n\pu). (A10) 



Using Eq. |A10[ we can calculate the expected value of 
Pn and its second moment, respectively denoted (Pn) 
and <P U ): 



(Pn 



PA) = 



Pn tt(Pn\n) dpn 



Pn Q{pn\n) dpn 



n + l 

(n + l)(n + 2) 



The standard deviation of Pn is therefore 



(P?i) - (Pn) 



Vn+1 



N, 



H 



Hence, given n measured threefold coincidences and 
Nh heralding signals, this conservative method yields 



J.CE) 

Pn 



n + 1 Vn+l 



N, 



H 



N, 



H 



where CE stands for "conservative estimation". In our 
results, we estimate the concurrence with both the max- 
imum likelihood and the conservative estimations of Pn. 



3. Experimental results 

The estimation of pooi Pio + Poi an d Pn proceeded as 
follows. For a time T — 166 hours, we let the interferom- 
eter drift to randomize the phase. During that period of 
time, we recorded the total number of heralded twofold 
coincidences between the heralding detector and detec- 
tor 1 (Nujj), and between the heralding detector and 
detector 2 (N 2 \h), and calculated N 2 = N^ H + N 2 \h- 
We also recorded the total number of heralded threefold 
coincidences, denoted N 12 \h, between all three detectors. 
We used a coincidence time window of 10 ns; this is the 
same value we used in the estimation of the concurrence 
based on the cross-correlation. We calculated 

N 2 y/N^ 



For pn, we used both the MLE, yielding 



(MLE) 

Pn 



2 27 | Nl2]H ± ^ Nl2 \ H 



N 



H 



N. 



H 



and the CE: 



JCE) 

Pn 



2.27 



N 12 \ H + 1 y/N 12 \ H + l 



N H 



N h 



Finally, we estimated poo = 1 — Pio — Poi ~ Pn- 

With a 10 ns coincidence window, we obtained p 10 + 
poi = 1.7777(34) x 10" 4 , V = 96.5(1.2)%, and observed 
two threefold coincidences (N 12 \h — 2). This yields 

p (MLE) = 2 g(2 1} x 1Q _9 

C (MLE) > 6.3(3.8) x 10- 5 

and 

p^E) = 3_g(2.2) x 10- 9 
C (CE) > 3.9(3.8) x 10~ 5 

Both the MLE and CE yield a lower bound for the con- 
currence that is greater than by at least one standard 
deviation. 
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